Quantum RLC circuits: charge discreteness and resonance 
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In a recent article [l| , we have advanced a semiclassical theory of quantum circuits with discrete 

charge and electrical resistance. In this work, we present a few elementary applications of this 

theory. For the zero resistance, inductive circuit, we obtain the Stark ladder energies in yet another 

00 I way; and generalize earlier results by Chandia et. al 3, for the circuit driven by a combination d.c. 

fy^ ■ plus a.c. electromotive force (emf). As a second application, we investigate the effect of electrical 

^^ ' resistance, together with charge discreteness, in the current amplitude, and resonance conditions 

Cn \ of a general RLC quantum circuit, including nonlinear effects up to third order on the external 

sinusoidal emf. 
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PACS numbers: 73.21.-b, 73.23.-b, 73.63.-b 



I. INTRODUCTION 



In a series of articles Li and Chen d, [j| and us [Ij B S B @ > have developed a theory of quantum electrical systems, 
^ ■ based on a treating such systems as quantum LC circuits; that is, electrical systems described by two fenomcnological 
parameters: an inductance L, and a capacitance C. Such quantum theory of circuits is expected to apply when 
the transport dimension becomes comparable with the charge carrier coherence length, taking into account both the 
quantum mechanical properties of the electron system, and also the discrete nature of electric charge. Now, in a 
recent work [1], we have proposed a semiclassical theory of quantum electrical circuits. This has been done with two 
G ■ goals in mind. First, to obtain useful predictions of the theory from very simple calculations, obtaining for example 
1-^ , energy spectra, and, second, to push the circuit analogy one step further, generalizing the equations of motion to 
r^ I include electrical resistance. As it is well known, the problem of electrical resistance, at the mesoscopic level, is that 
Q ' of the contact resistance, i.e., the coupling between the mesoscopic system and a macroscopic system with which it 
O . interacts, which leads to the conductance quantization condition, as it was first shown by Landauer [lO|. 

The semiclassical theory of quantum RLC circuits [l| starts from the quantum Hamiltonian of the LC circuit [l|, 
T— I ' SH, BB 0, HQi ^'^'i adds a term to account for electrical resistance, R(j)/L. The resulting equations become 

>'. 

rn ' q = -^sin(0/(/)o). 

• ■ L/ 

!>■ 

f^ ' The equations above are considered, mathematically, as classical equations, but they include quantum effects, the 

00 , quantized nature of electric charge (through the parameter 0o = fi-/Qe), and electrical resistance R; which couples the 

^^ ' system to the external reservoir Jj]. These equations are highly nonlinear, but they reduce to the usual equations of 

J> I the RLC circuit, in the discrete charge limit, q^ -^ 0. An interesting consequence is that it is possible to find stable, 

constant charge and flux {cf) = 2mT(j)o) solutions of the equations [I] such that, if one further assumes that the charge 

is discrete {q — —mq,,), then the following relation should exist between the system parameters jlj, which bears a 

S I strong resemblance to the Landauer formula [lO| 

_ RC _ mql_ _ m 
L n h n 

In this work, we apply our theory to some simple systems; in section |TT1 we consider the so-called L -design model, a 
quantum circuit under a constant electromotive force e, obtaining the Stark ladders in a very simple manner [2|. Next, 
in section IIIIl we consider the same L design case, add electrical resistance, according to our semiclassical theory. 



•Electronic address; lcutreras@uach.cll 



and subject it to a combined a.c. plus d.c. electromotive force, e{t) = Eq + Si cos(wi). In section ITVl we consider a 
quantum RLC circuit under a sinusoidal-only external electromotive force e(t) = £oCOs(tL'i), and study the resonance 
conditions, considering nonlinear effects up to £q, and electrical resistance. We find that the resonance frequency 
of the quantum dissipative system is lowered, with respect to the non dissipative system, the shift is found to be 
Ao^/wo = —{£o/R)'^/{4:(f>Q/L)'^. We remark that the results obtained show the usefulness of the semiclassical approach 
to both dissipative and non-dissipative quantum circuits, therefore, it provides a way to advance our understanding 
of quantum circuits, and encourages us to pursue the subject in future works. 

II. INDUCTIVE CIRCUIT: STARK LADDERS 

In a recent work, Chandia et. al. [4I, using the quantum circuit approach have obtained the energies of the Stark 
ladders, previously conjectured by Chen at. al. [ll|. Chandia et. al. [3| consider a quantum circuit without resistance, 
in the L design situation (the capacity, C -^ cxd) situation, under a constant applied electromotive force e. In this 
work, we that the semiclassical method described in a previous article [l| may be used to obtain the quantized 
energies of the stark ladders. To prove this, consider a quantum circuit (and large, C ^ 00 capacity), under a d.c. 
electromotive force e. Using 0o — ^/<le, the equations of motion may be written as 

I = ^sHHM, (4) 

and the Hamiltonian is 

H=-qe+-^sm\(f>/2c^o)- (5) 

To find the energy spectrum, we use the semiclassical method [l| , finding first the charge q as a function of the energy, 

g=^(-i?+^sin2(<^/20o)). (6) 

Next, we compute the action variable J, 

. = /,.* = J^^-%.,.-?^(^-f), (7) 

and impose the quantization condition in the semiclassical way (Bohr-Sommerfeld), J — nh, we obtain precisely the 
result of Chen et al [111], also obtained by Chandia et al. [^, namely 

E=-nqee + -^. (8) 

It it interesting to see that these results are very robust, in the sense that they are preserved, under very different 
mathematical descriptions, depending strongly on the charge quantization condition more that anything else. 

III. INDUCTIVE CIRCUIT: COMBINED A.C, D.C. AND RESISTANCE 

In their recent article, Chandia et. al [2| consider a quantum inductive, nondisipative circuit, subject to a weak 
external electromotive force given by e{t) = Eq + Ei cos(cji). They showed that the (physical) electrical current in such 
a quantum circuit, in the first order approximation, has a nonzero time average over a period T = 2tt/ujb, given by 

dq EiU)B ( ii\v?{-KUj/ujB) \ . . 

In this equation lob = QeEo/h = Eq/(J)q (0o = J^/Qc is a flux quantum). To arrive to equation ([9|), one computes the 
flux, replaces into equation Q, and expands the sine function up to first order, assuming that the alternating voltage 
£1 is small, i.e. Ei << u!4>o, 



dq (f>a f . , ,s , SiujB . / ,N / ,N 

—r — -r \ smfwet) H smlujt) cosiojBt) 

at -L \ Sqi^ 



(10) 



Equation ^ is obtained after averaging equation ©over a time period T = 27^ /ujb- The results !2J|, tell us that 

• The electric current has zero average when uj = ulob, for integer n. 

• The electric current has nonzero average for ui ^ nuiB, and it shows extrema (resonances) when oj « (n + l/2)ujB; 
the strongest maximum occurs for lu rs lobI'2- 

• For small frequencies, lo « lub, 
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(11) 



Another important observation is that, in alternating current circuits, it is customary to average with respect to the 
pariod of the source, and also that the observable quantity would be the root mean square of the current, ■\/< (p >, 
instead of the average value of the current. First, define the parameter k = eiUJB/i^o^)^ then compute the average, 
for the square of the current, given below, neglecting k^ terms, 

9^ = (— )^ (sin^iLUBt) + ksm{2ujBt)sm{ujt)) . 



The results are. 
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• Average for T = 27r/a; 
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(13) 



Now, let us study the changes that may be observed when one includes electrical resistance in our description of the 
inductive, discrete-charge, quantum circuit discussed previously. As it has been shown by previously ^J, the equations 
that describe an inductive circuit with electrical resistance under the emf eo + £i cos(wt) are 



it , , 

+ — = eo+eicos{ujt) 

q = sm( — ). 



(14) 
(15) 



Since these equations remain uncoupled, they may reduced to an integration. The particular solution for cf) = 4>p{t) 
may be written as 



4>p(t) = A + B cos{ujt - a). 



(16) 



Use the complex impedance Z{uj) = R + iojL, and |^(w)| — \/B? + {loLY, and the phase angle a, tan(a) = ujL/ R, 
then the coefRcients A and B become 



A=^ (17) 

The solution obtained so far docs not take into account the initial condition at t = 0, therefore we must add an 
exponentially decaying term, e"*'"^, r = L/R, so that 

(^(t) =yl + Bcos(wi-a) + Ce~*/^. (19) 

The constant C is chosen so that (j){t = 0) = 0, then C = —A — B cos(a), then the proper solution is 

0(t) ^ a(i- e^*/^] + B {cos{ujt ~a)- cos(a)e-*/^] . (20) 

This solution has the correct behaviour as i? ^ 0, namely, it coincides with the solution from the R = case, equation 
(|3]); however, there are significant differences, since: 

• The solution (j){t) for the case R — increases linearly with time, while, for R ^ 0, (pit) grows linearly with time 
only for t << t 



• For t >> T, and R =/= 0, the solution (j){t) is always bound, and it coincides with ppit), in equation (jl6[) . 

We conclude that the effect described by Chandia et. al [2| may only be observed for very weakly dissipative systems, 
for short times, t << t, so that the observation times involved {^-k Jlo and 2ti/lob) should be much shorter than the 
characteristic time r = L/R. For longer times, t » t, the current will be 

00 . f Leo Lei , , \ ,„,^ 

q=— sm — - + cos wi - ") ■ 21) 

L \R(pQ \Z{u})\(\)Q J 

Now, let us consider the case in which Lei/\Z{ijj)\(f)Q « 1, then 

/ T,Pr. \ 1 

(22) 
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so that the time average of the electric current, over a period of the source, T — 2Ti/bj, is 

Notice that the current q has nonzero time average, except for some particular values of the argument of the sine 
function, Leo/Rfjjo = mr, in other words, Eq — wk^q/t, a very simple result. 

Now, for i >> T, the solution for q is given by equation (P^ . therefore, averaging over a period of the source 

(T = 27r/w) 

IV. RLC CIRCUIT: RESONANCE 

Consider now the circuit under a purely sinusoidal perturbation, e{t) = e^ cos{ujt), and now look for a solution as a 
pertubative series, we obtain only the first two terms of the series. We write down the circuit equations once again. 



R , q , . 

— 0+ — = eocos(ujt) 



90 ■ 
— sm 



Let us seek solutions as 



q{t) = qi{t) + qs{t) + ■ ■ ■ 

we write them like this since we know that upon series expansion the solution qi (and 0i) will be linear on the 
parameter Eq, and the solution ^3 (and (p^) will be proportional to Eq, with only odd-order terms appearing. Now, 
insert the series solutions into the RLC equations, obtaining the linearized equations 

(25) 

(26) 

(27) 

(28) 

It is easy to see that both equations may be solved by the complex-number method. Therefore it is convenient to 
define the complex impedance Z{uj) — R + iuL + l/itoC = |Z(w)|e*"% as well as its phase angle and magnitude, 

Z{u;) = R + iujL + l/iujC (29) 

\Z{ij)\ = ^R^ + {ujL - l/wC)2 (30) 

LoL — 1/ujC 

tan(ai) = -^ . (31) 

K 

In this way, the first order pseudo fiux (pi and pseudo current are 

I,it) ^ M!) ^ ^^ cos(^t - «i)- (32) 

L \Z{uj)\ 





= eocos(tjt) 
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To write the current q, we define the parameters, 

(0?)^ L' el 



HlL 60§ \Z{uj)\^ 



(33) 



tan(a3) = -^ (34) 

tan(/3i) = ^ (35) 

tan(/33) = ^ (36) 

The current q below contains terms oscilating at the frequency w, and also at Sw, which comes from the cubic 
nonlinearity. 



«(*) = \7( \\ cos{ujt-ai) '' cos{ujt ~ 2ai + Pi) + 

\Z{uj)\ MZHl 

4\Z(3uj)\ cos(3wi - 3ai - a3 + /^s). 



The results of the previous paragraphs allow us to see that there is a slight change of the resonace frequency for 
a quantum circuit, the corrections are small, and depend in a nonlinear way on the external source amplitude eg. 
Define the dimensionless parameter k << 1, 



im^R-^ + {loLY _ v/i?2 + {ujLfiLeof 



4en 



mzi^w 



(37) 



the terms oscilating at frequency lo may be put together in a single term, then the current may be written as (defining 
the additional phase angle A0). We neglect the term at frequency 3w, since if one evaluates its root mean square 
average for a period of time T = 2t: /uj, it will turn out a zero value. 
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So 



\Jl + K? — 2k cos(— ai + j3x) co8{u}t — ai + /S.9). 
[1 — KCOs(— ai + /3i)] cos(a;i — ai + A0) 



/o(u;)[l - A{ijj)] cos{ujt - ai + A6') 



where we have defined the constants 



loH 

A{oj) 
I{u) 



£0 



\Z{u)\ 

KC0s(/3i — ai) 
h{l-A{uj)). 



(38) 

(39) 
(40) 

(41) 

(42) 
(43) 



We may put out main result, /(w), in simpler terms, to do that, use the addition formula for cosine, and the definitions 
given earlier for the paramenters j3i, a\, etc.; after some algebra, we get 



cos(/3i — ai) 
A{l,) 



R^ +ujL(ujL- \/ujC) 
{Leo? 



Hi\Z{Lo)\'^ 



R^ +ujL{ljL- I/cjC)] . 



(44) 
(45) 



From this we shall compute the shift on the resonance frequency due to both charge quantization and electrical 
resistance. To do this, we expand the equation above about lj ^ luq ~ 1/^/LC, up to second order. 



A{lo) 



(Leo? 



blR^ 



2L{lu-luo) 7L'^{uj-u)o? 



/o(-) = I 



1- 



LuoCR^ 
2L^iLU - Up? 

i?2 
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(46) 
(47) 



We need to solve dl/duj — 0, for Aw — oj — ujq, recalling that A{uj) « 1, and that dIo{uJo)/doJ — 0, 



dl{u 



duo 
dl(uj) 



dlo(uj) ._ ,, ^. -, , .dA(uj) 



duo 



dio duo"^ 

using the expansions above, we obtain the solution 



dio 



Auj 



\ duo 



{eo/Rfm 



duj^ 



IHIR'' I6(0o/i)^ 



(48) 
(49) 

(50) 



V. FINAL REMARKS 

In this article, we have studied some simple applications of the ideas of discrete-charge quantum RLC circuits, and 
applied to two simple systems. In the first place, we reobtained the energy spectrum (Stark ladders) for the model 
L-design system, using the semiclassical method [i|. Next, we consider the L design case studied by Chandia [4], 
but including the effect of external resistance, finding that their results hold under more stringent conditions, namely 
that the observation times should be much shorter that the characteristic time r = L/R] finding the correction in 
the other case. Finally, we used the generalized quantum RLC circuit equation to study the RLC circuit under a 
sinosoidal electromotive force. We have shown that the resonance frequency of the system becomes shifted, and the 
shift is proportional to the square of the ratio of a classical current Eq/R and a quamtum current (t>o/L; a simple result 
may be tested under appropiate experimental conditions. 
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